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Abstract 



We investigate the eigenvalue spectrum of the staggered Dirac matrix in SU(3) 
gauge theory and in full QCD on a 6 3 x 4 lattice. As a measure of the fluctuation 
properties of the eigenvalues, we study the nearest-neighbor spacing distribution 
P(s) for various values of (3 both in the confinement and in the deconfinement phase. 
In both phases except far into the deconfinement region, the lattice data agree with 
the Wigner surmise of random- matrix theory which is indicative of quantum chaos. 
We do not find signs of a transition to Poisson regularity at the deconfinement phase 
transition. 
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1 Introduction 

> 

•rH , 

The properties of the eigenvalues of the Dirac operator are of great importance 
for the understanding of certain features of QCD. For example, the accumu- 
lation of small eigenvalues is, via the Banks-Casher formula [1], related to the 
spontaneous breaking of chiral symmetry. Recently, the fluctuation properties 
of the eigenvalues in the bulk of the spectrum have also attracted attention. It 
was shown in Ref. [2] that on the scale of the mean level spacing they are de- 
scribed by random-matrix theory (RMT). In particular, the nearest-neighbor 
spacing distribution P(s), i.e., the distribution of spacings s between adja- 
cent eigenvalues on the unfolded scale (see below), agrees with the Wigner 
surmise of RMT. According to the so-called Bohigas-conjecture [3], quantum 
systems whose classical counterparts are chaotic have a nearest-neighbor spac- 
ing distribution given by RMT whereas systems whose classical counterparts 
are integrable obey a Poisson distribution, P(s) = e~ s . Therefore, the spe- 
cific form of P(s) is often taken as a criterion for the presence or absence of 
"quantum chaos". It should be noted, however, that there is no universally 
accepted proof of the Bohigas-conjecture yet. The field of quantum chaos is 
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still developing, and there are many open theoretical problems. For reviews 
on this subject we refer to Ref. [4]. 

A transition in P(s) from Wigner to Poisson behavior is observed at the 
localization transition in disordered mesoscopic systems [5]. In Ref. [2], the 
question was raised if there is such a transition from chaotic to regular behavior 
in the case of the lattice Dirac operator, and the present study serves to 
investigate this question. For a preliminary account of results for quenched 
QCD, see Ref. [6]. Recently, a Wigner to Poisson transition was also studied 
in the context of a spatially homogeneous Yang-Mills-Higgs system [7]. The 
question if chaos in an iV-component <^ 4 -theory in the presence of an external 
field survives quantization was investigated in Ref. [8] . An analytical treatment 
of the transition between Wigner and Poisson behavior can be found in Ref. [9]. 

In RMT, one has to distinguish between several universality classes which are 
determined by the symmetries of the system. For the case of the QCD Dirac 
operator, this classification was done in Ref. [10]. Depending on the number of 
colors and the representation of the quarks, the Dirac operator is described by 
one of the three chiral ensembles of RMT. As far as the fluctuation properties 
in the bulk of the spectrum are concerned, the predictions of the chiral ensem- 
bles are identical to those of the ordinary ensembles [11]. In Ref. [2], the Dirac 
matrix was studied in SU(2) using both staggered and Wilson fermions which 
correspond to the symplectic and orthogonal ensemble, respectively. Here, we 
study SU(3) with staggered fermions which corresponds to the chiral unitary 
ensemble. We thus cover the last remaining symmetry class in the confined 
phase. Although the unitary ensemble is, from a mathematical point of view, 
the simplest one, the RMT result for the nearest-neighbor spacing distribution 
is still rather complicated. It can be expressed in terms of so-called prolate 
spheroidal functions, see Ref. [12] where P(s) has also been tabulated. A very 
good approximation to P(s) is provided by 

P(s) = . (1) 

This is the Wigner surmise for the unitary ensemble. 



2 Eigenvalue analysis 

We generated gauge field configurations using the standard Wilson plaquette 
action for SU(3) with and without dynamical fermions in the Kogut-Susskind 
prescription. We have worked on a 6 3 x 4 lattice with various values of the 
inverse gauge coupling (3 = Q/g 2 both in the confinement and deconfinement 
phase. The boundary conditions were periodic for the gluons and periodic in 
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space and anti-periodic in Euclidean time for the fermions. We typically pro- 
duced 10 independent equilibrium configurations (separated by 1000 sweeps) 
for each /3. Because of the spectral ergodicity property of RMT one can replace 
ensemble averages by spectral averages if one is only interested in bulk proper- 
ties, and thus one does not need a large number of independent configurations 
to compute P(s). 

By definition, the Dirac operator p = <f) + ig^ from which the spectrum 
is calculated does not contain the quark mass term. The operator is anti- 
hermitian so that all eigenvalues are imaginary. For convenience, we denote 
them by iX n and refer to the X n as the eigenvalues in the following. Because 
of {^, 75 } = 0the X n occur in pairs of opposite sign. All spectra were checked 
against the analytical sum rules 

£A n = and £A2=3^, (2) 

n \ n >0 

where V is the lattice volume. 

To construct the nearest-neighbor spacing distribution P(s) from the eigenval- 
ues, one first has to "unfold" the spectra. This procedure is a local rescaling 
of the energy scale so that the mean level spacing is equal to unity on the 
unfolded scale. One first defines the staircase function N(E) to be the number 
of eigenvalues with A < E. This staircase function has an average part and 
a fluctuating part, N(E) = N aY (E) + N R (E). The average part is extracted 
by fitting N(E) to a smooth curve, e.g., to a low-order polynomial. One then 
defines the unfolded energies to be x n = N av (E n ). As a consequence, the se- 
quence {x n } has mean level spacing equal to unity. Ensemble and spectral 
averages are only meaningful after unfolding. 



3 Results and Discussion 

We begin the discussion of our results with the pure gluonic case. Figure 1 
shows the staircase function N(E) which in the continuum is given by N(E) = 
Jq p(X)d\, where p(A) is the spectral density. Note that p(A) cannot be ob- 
tained universally from a random-matrix model. The staircase function is not 
very sensitive to the value of /3, neither in the confinement nor in the decon- 
finement regime. However, across the phase transition N(E) is diminished for 
small E reflecting the suppression of small eigenvalues. 

Figure 2 shows the nearest-neighbor spacing distribution P(s) corresponding 
to the parameters of Fig. 1, compared with the RMT prediction. In the confine- 
ment phase, we find the expected agreement of P(s) with the Wigner surmise 
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Confinement 



Deconfinement 



(5 = 2.8 (solid line) 
P = 5.0 (dashed line) 




(5 = 6.0 (solid line) 
P = 10.0 (dashed line) 




Fig. 1. The staircase function N(E) representing the number of positive eigenvalues 
< E for typical configurations on a 6 3 x 4 lattice in pure SU(3). 



of Eq. (1). In the deconfinement phase, we still observe agreement with the 
RMT result up to (3 = 10.0 (we have also plotted the Poisson distribution for 
comparison) . 

We continue with the case of full QCD with Nf — 3 degenerate flavors of 
staggered quarks with mass ma = 0.1 and ma = 0.05, respectively. The stair- 
case function N(E) in Fig. 3 shows practically no dependence on the masses 
considered. 



Figure 4 compares the nearest-neighbor spacing distribution P(s) of full QCD 
with the RMT result. No mass dependence is visible. Again, in the confinement 
as well as in the deconfinement phase we still observe agreement with the RMT 
result up to (3 = 8.0 (not shown). The observation that P(s) is not influenced 
by the presence of dynamical quarks could have been expected from the results 
of Ref. [11]. Those calculations, however, only apply to the case of massless 
dynamical quarks. Our results, as those of Ref. [2], strongly indicate that 
massive dynamical quarks do not affect P(s) either. 

No signs for a transition to Poisson regularity are found. Thus, the decon- 
finement phase transition does not seem to coincide with a transition in the 
spacing distribution. For very large values of f3 far into the deconfinement 
region (not shown), the eigenvalues start to approach the degenerate eigen- 
values of the free theory, given by A 2 = X^=i si n2 (27rn^/L^)/a 2 , where a is 
the lattice constant, is the number of lattice sites in the /z-direction, and 
n M = 0, . . . , L M — 1. In this case, the nearest-neighbor spacing distribution is 
neither Wigner nor Poisson. However, it is possible to lift the degeneracies of 
the free eigenvalues using an asymmetric lattice where L x , L y , etc. are rela- 
tive primes [13]. For large lattices, the nearest-neighbor spacing distribution 
of the non-degenerate free eigenvalues is then given by the Poisson distribu- 
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Confinement Deconfinement 

£ = 2.8 £ = 6.0 




Fig. 2. The nearest-neighbor spacing distribution P(s) averaged over 10 independent 
configurations (except for (5 = 2.8 where all single configurations are plotted to give 
an indication of the deviations) on a 6 3 x 4 lattice in pure SU(3) (histograms) 
compared with the random-matrix result (solid lines). For comparison, the Poisson 
distribution P(s) = e~ s is plotted for f3 = 10.0 (dotted line). There are no changes 
in P{s) across the deconfinement phase transition. 



Confinement (3 = 5.2 



Deconfinement (5 = 5.4 



ma = 0.1 (solid line) 
ma = 0.05 (dashed line) 




ma = 0.1 (solid line) 
ma = 0.05 (dashed line) 




Fig. 3. The staircase function N(E) for typical configurations on a 6 3 x 4 lattice in 
full QCD with different quark masses. 
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Confinement j3 = 5.2 



Deconfinement j3 = 5.4 



ma = 0.1 ma = 0.1 




Fig. 4. The nearest-neighbor spacing distribution -P(s) averaged over 10 independent 
configurations on a 6 3 x 4 lattice in full QCD (histograms) compared with the 
random-matrix result (solid lines). Again, there are no changes in P(s) across the 
deconfinement phase transition. 

tion. While it may be interesting to search for a Wigner to Poisson transition 
on such asymmetric lattices, it is unlikely that such a transition will coincide 
with the deconfinement phase transition. 

We do not believe that the absence of a signature for a transition from Wigner 
to Poisson behavior at the deconfinement phase transition is due to the finite 
lattice size. Even for the small lattice size we used, the agreement of P(s) with 
the RMT curve is nearly perfect. This leads us to believe that we should have 
seen some sign of a transition if it existed in the thermodynamic limit. 



4 Conclusions 

We have searched for a transition in the nearest-neighbor spacing distribu- 
tion P(s) from Wigner to Poisson behavior across the deconfinement phase 
transition of pure gluonic and of full QCD. Such a transition exists, e.g., at 
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the localization transition in disordered mesoscopic systems [5]. In a Yang- 
Mills-Higgs system a smooth transition along a Brody distribution was seen 
[7]. We observed no signature of a transition in our data for the lattice Dirac 
matrix, neither for pure SU(3) nor for full QCD. The data agree with the 
RMT result in both the confinement and the deconfinement phase except for 
extremely large values of f3 where the eigenvalues are known analytically. Our 
analysis of full QCD with two different quark masses showed no influence on 
the nearest-neighbor spacing distribution. 

Our findings are consistent with earlier studies of the chaotic dynamics in 
classical lattice gauge theories which have shown that lattice gauge theories 
are chaotic as classical Hamiltonian dynamical systems [14]. Also, it was found 
recently that the leading Lyapunov exponents of SU(2) Yang-Mills field con- 
figurations indicate that configurations in the deconfinement phase are still 
chaotic, although less chaotic than in the strong coupling phase at finite tem- 
perature [15]. 

In hindsight, our results are not totally unexpected. Temporal monopole cur- 
rents survive the deconfinement phase transition leading to confinement of 
spatial Wilson loops. Thus, even in the deconfinement phase, the gauge fields 
retain a certain degree of randomness. It was shown that the relevant parame- 
ter that drives the transition from Poisson regularity to chaos is of the order of 
the mean level spacing A of the system [16,9]. One could search for a formal 
separation of the Dirac operator into a regular and a chaotic contribution, 
p =ftrcg + (oc/ A) -^chaotic- Then, a chaotic perturbation of size a~ A will gen- 
erate a Wigner distribution. In future investigations, it would be interesting 
to try to disentangle the corresponding contributions to the Dirac matrix. 
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